clc;
clear;
f=@(x)exp(-x.^2);
I=integral(f,0,inf,'RelTol',1e-20)
syms x;
I1=int(exp(-x^2),0,inf)
vpa(I1)